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摘要 : 气 - 液 两 相 驱 替 在 油气 藏 开发 中 广泛 存在 。 传 统 的 渗流 理论 认为 压力 和 达 西 速度 在 非 连续 宽 面 处 是 连续 
的 , 我 们 在 以 往 的 研究 对 其 进行 了 证 伪 。 本 文 在 前 期 不 可 压缩 流体 的 基础 上 ， 以 气 -水 两 相 驱 蔡 为 例 , 研究 了 多 孔 介 
质 中 移动 的 气 -水 突变 界面 处 的 压力 和 渗流 速度 的 连接 条 件 。 研 究 结果 表明 ， 在 气 -水 两 相 驱 蔡 界 面 处 : (1) 总 体 压 
连续 ， 但 流体 压力 可 以 不 连续 ;(2) 总 体 流速 的 分 布 可 以 不 连续 ， 但 满足 特定 的 函数 关系 ;“3) 流体 各 相 渗流 速度 
不 连续 。 

Abstract: Gas-liquid two-phase displacement are widely observed in oil and gas reservoir. The traditional 
seepage theory holds that the phase pressures and Darcy velocities must be continuous at the jump interface, 
and we have falsified them in previous studies. In this paper, the jump interface condition of two-phase seepage 
is extended from incompressible fluids to compressible. Taking gas-water two-phase displacement as an example, 
the new connection conditions of phase pressures and Darcy velocities at the gas-water displacement front in 
> porous media were built. The results show that at the flood front, the global pressure is continuous, but not the 
(© fluid pressure are. (2) The total Darcy velocity can be discontinuous with a specific functional relationship;(3) The 
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© phase Darcy velocity of each phase of the fluid is discontinuous. 
S 
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= 工 介绍 
D TRA AE SS CET UE AC Piz FEE, 例如 : 边 底 水 类 型 气 藏 、 含 水 层 型 地 下 储 气 库 库 气 库 、 
SS CO2 地 下 埋 存 、 注 气 驱 油 、 注 气 驱 油 、 带 气 顶 油 茂 开发 天然气 排 水 成 茂 、 动 态 法 测试 气 -水 两 相 滩 
Se 透 率 等 储 多 问题 。 气 - 液 两 相 驱 普 的 连接 条 件 是 建立 渗流 数学 模型 ， 跟 踪 界 面 移动 ， 计 算 驱 普 作 用 下 
O 产量 。 根 据 渗流 文献 ， 在 驱 蔡 前 缘 界 面 处 流体 各 相 压力 和 滩 流 速度 连续 。 即 ; 
二 V, = Vp V, =V, 
© P= 0-2, (1) 
P. =P. 


根据 质量 守恒 方程 , 单位 时 间 内 微 元 内 流体 变化 F 流 入 微 元 的 流量 | 流出 微 元 的 流量 |。 对 于 移动 界 


面 ， 体 积 为 0， 没有 源 汇 项 ， 则 有 ， 了 流入 微 元 的 流量 | 流出 微 元 的 流量 F0。Hassanizadeh, S. M., & Gray, 


W. G. (1989) 用 物质 守恒 方程 对 其 进行 了 证 明 。 针 对 界面 两 侧 区 域 毛管 力 不 同 的 情况 ， 作 者 在 文献 中 
对 传统 的 界面 条 件 进行 质疑 , 然后 根据 物质 守恒 方程 和 两 相 渗流 公式 推导 了 新 的 界面 条 件 以 及 路 突变 
界面 的 渗流 公式 中 。 气 - 液 两 相 驱 蔡 界 面具 有 三 个 特征 : (1) 驱 蔡 界面 随时 间 不 断 移 动 ; (2) 驱 蔡 界 
面 一 般 是 突变 界面 ， 流 体 分布 可 能 不 连续 ; G) 气体 压缩 性 高 ， 驱 替 过 程 需要 考虑 压缩 性 。 所 以 本 文 
将 在 前 期 不 可 压缩 流体 两 相 突 变 界 面条 件 的 研究 基础 上 ,考虑 流体 的 压缩 性 , 研究 多 孔 介质 中 不 互 溶 
气 - 液 两 相 驱 蔡 前 缘 的 界面 条 件 。 


2 质疑 对 传统 突变 界面 条 件 主要 论点 
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(1) 驱 蔡 前 缘 界 面 的 界面 条 件 需要 建立 的 侧面 两 侧 不 同 流体 质点 的 压力 连续 条 件 和 流体 速度 的 
连接 条 件 如 图 ); 体积 为 0 的 医 面 微 元 流入 吕 流 册 量 所 指 对 象 是 在 界面 上 同一 质 车 说 明 两 者 描述 的 
不 是 同一 物理 过 程 ， 


: : 
二 + 
X Xp x 


1 界面 条 件 与 界面 上 质点 守恒 方程 所 指 对 象 


(2) 以 一 维 为 例 ， 假 如 包括 界面 的 薄 层 状 微 元 厚度 为 b>0， 流 入 微 元 的 流速 记 为 ww ， 流 出 微 元 
的 流速 记 为 vi (= o,w)， 两 者 的 数量 关系 即 为 连接 条 件 , 应 满足 根据 质量 守恒 方程 (2), 根据 该 方程 并 
不 能 得 出 v; =v: 而 传统 的 渗流 力学 针对 界面 条 件 所 用 的 质量 守恒 方程 表达 为 (3) ” ,方程 (2) 与 (3) 
的 区 别 实质 上 已 经 蔡 换 了 质量 守恒 方程 。Hassanizadeh, S. M., & Gray 试图 证 明 由 (2) 可 以 得 到 (3), 
但 我 们 发 现 证 明 过 程 有 误 。 


0 0 
ap Pelt a 0)=0 How) r 
paki A (1=0,w) 
dx 
py i (¢ =O, w) (3) 
Py =P: 


2 包括 驱 蔡 界面 的 江 薄 层 状 微 元 示意 图 
(3) Buckley-Leverett 模型 的 解析 的 解 明确 地 显示 渗流 速度 在 前 缘 突 变 , 与 公式 (2) 明 显 区 别 ; De 
Neef, M. J., & Molenaar, J. (1997) 公 式 〈1) 显示 流体 各 相 压 力 和 毛管 力 不 能 同时 满足 ， 为 此 提出 改进 
的 界面 条 件 ， 他 们 认为 润 相 流 体 压力 连续 ， 但 非 润 湿 相 压力 可 以 不 连续 。 
(4) 对 于 移动 界面 而 言 , 公式 (1) 所 示 的 传统 的 界面 条 件 与 Rankine-Hugoniot 条 件 矛 盾 ; Rankine- 
Hugoniot 条 件 可 由 质量 守恒 方程 严格 推导 。 以 公式 (2) 为 例 ， 其 Rankine-Hugoniot 条 件 为 : 
v -=v v =v} 
“A-s Alsi) a 
(5) 因为 按照 公式 (1) 所 示 的 传统 的 界面 条 件 ， 以 及 两 相 流 达 西 公式 ， 突 变 界面 处 传导 率 应 该 为 
界面 两 侧 的 调和 平均 值 ， 它 不 能 正确 反映 真实 的 渗流 过 程 ， 所 以 被 单 点 上 游 权 〈SPU) 和 其 它 高 阶 上 
游 权 方法 取代 。 调 研 文件 我 们 发 现 虽然 (1) 表 示 的 界面 条 件 广泛 见 诸 文 献 ， 但 并 没有 得 到 真正 应 用 。 
27 7 ee, ee 
Ip +I M Ax m Ax 
(6) 一 些 学 者 能 通过 和 否定 流体 饱和 度 不 可 能 呈 间 断 分 布 ， 从 而 解释 突 主 界面 处 压力 和 流速 必须 


(pr -pi),T = (£=w,g) (5) 


连续 。 我 们 的 观点 认为 : 首先 从 微观 而 言 ， 流 体 分 布 界面 可 以 不 连续 ， 宏 观 是 微观 的 统计 ， 所 以 在 宏 
观 上 流体 的 饱和 度 能 呈 突 变 分 布 ， 其 次 ， 从 工业 计算 角度 ， 突 变 界面 区 域 对 于 整体 的 渗流 阻力 很 小 ， 
压力 降 主 要 发 生 于 突变 界面 两 侧 ， 即 使 流体 分 布 以 连续 分 布 的 方式 突变 ， 通 过 转换 为 不 连续 分 布 ， 可 
得 到 相近 的 结果 ， 而 且 可 大 幅度 减 小 计算 量 。 


0 iti 3 
图 3 将 流体 处 突变 但 连续 分 布 界面 简化 为 不 连续 界面 示意 图 
Ce: 真实 界面 附近 流体 分 布 ， 右 : 简化 后 的 不 连续 界面 ) 


3 不 互 溶 气 -水 两 相 渗流 方程 
以 不 互 溶 的 气 -水 两 相 为 例 根据 气 藏 渗流 力学 ， 有 气 -水 两 相 渗 流 的 质量 守恒 方程 : 
ð 


b— (PS) +V (P) =0 
0 
Pa (Pr)tV (prs)=0 
S, +5, =1 
ia (6) 
V | | 


deele 
ôt| B, A 


Pe 5 P7 Pw (7) 


Bs 表示 地 层 体 积 系数 ， 无 因 次 ; ZANT EAT TAR: S% 表 示 含 水 饱和 度 ; v 表示 达 西 
流速 m/day. OFAN x 的 水 平 夹 角 。 流 体 在 气 藏 中 附 合 达 西 定律 ， 所 以 有 : 


v, = (Vp,—p,g sind) (l=w,g) (8) 
4 气 - 水 两 相 驱 蔡 前 缘 的 界面 条 件 


4.1 驱 痊 前 缘 处 总 体 流速 的 连接 条 件 
以 一 维 坐标 为 例 ， 方 程 (7) 改写 为 


的。 
ôt| B, cr B, 
OEE, (5) 
ot\ B, Ox\ B, 


(9) 


考查 驱 蔡 前 缘 附 近 [x x, KR, 用 大 表示 气 -水 的 驱动 前 缘 界面 , 随时 间 移动 ，x >a, xe 


F + = + 。 
We, E, E, > 0,6, >0 BE, << 克 一 让 ,2 << r -r Hp ie =e rn =r+eé,, 即 有 : 


[7.-é,.7+6,1C[7,.7 | (10) 


驱 蔡 过 程 中 ,在 极 短 的 时 间 可 认为 @@ 驱 蔡 前 缘 界 面 以 形态 保持 不 变 的 方式 移动 ; @ 当 孔隙 度 连 续 分 布 ， 
可 认为 在 此 区 间 孔 际 度 为 常数 。 
对 方程 公式 (9) - (a) 两 端 在 [x.,] 范 围 内 积分 有 : 
f 8), dr+| | (1) 
» ôt| B, % Br B, 
Xp S Xp S xp S Xp S i 
“j g= dxf g= x+ | p-e x+ f g= |dx+|- =0 
Ot “xa B, Ot °F B, Ot “3 B, Ot 7È B, B, - 
(12) 
a Sr ae (S ae | .8. \a a 
COREG a 过 Gane w]e] =0 (135 
Se 和 ,为 常数 于 是 有 : 
Se s? + > 
| :| a | Aah -0 a4) 
B; Ot B; Ot B, be 
a b 
Ve _ Ve 
Be B° 
i: gs (15) 
Ot a d 
B, B, 
等 式 C15) 正 是 一 维 气 -水 两 相 驱 蔡 前 缘 界面 的 Rankine-Hugoniot 条 件 ， 渗 流 方程 的 ， 同 理由 〈9) 
- (b) 可 得 ; 
re 
a b 
Or Ba By) (16) 


a b 
ve. V 4 b 
8 8 Vw = Vi, 
Bo B° Bet Be 


8 8 = ŻW w (17) 


如 果 采 用 径 向 坐标 系统 仍然 可 以 得 到 等 式 (15)〉 ~ (17) 
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Al 4 边 底 水 径 向 驱 替 示意 图 


Ve Kk 


ENX: f ==- Vp,- p,gsinf) (l=w, (18) 
Ve B, Bus | Pp Ps ) ( W, 8) 
3 Ny a k v 
eM: A =— ,VD =, +0, ,=— (f=w, (19) 
| LB, om 5 me) 
4A 总 7 dP, x dp = 
ma ¥, =-KA,| —*-f,—* +7, (20) 
dx dx 
o,=f,0, +GP 
-T asl, Qn) 
V, = f,Y, —GP, 


e z x ; ; ~ | dp. 
其 中 : Fo = (fpsg + fps8)sing, Ay, =(p, -ps )8 sind 5 OP. = KA, (Zerav) 


将 (21) AR (17) 可 得 到 : 


a b 
vi=nv +B (22) 
其 中 : 
hee GP: GP“ GP? GP’ 
b b a B° B? B? 
esi- S? Ss- S? ge ay @——w ge w Se —S° g g1 gt 
BY B* B 8 B 8 B? B BY 8 g 
7 = = 
fo Fen 
a BS b a By b a B; b a By b 
Sa pe Se Swope Ss Sog Sw pe Sw 


(23) 


说 明 在 驱 蔡 前 缘 两 侧 ,， 总体 渗流 速度 不 相等 。 对 比 不 可 压缩 流体 可 知 在 驱 蔡 前 缘 两 侧 有 总 体 速 度 
相等 ， 说 明 不 可 压缩 流体 和 可 压缩 流体 具有 不 同 的 界面 连续 条 件 ， 前 者 也 可 以 看 作 是 后 者 的 特例 。 


4.2 驱 替 前 缘 界面 处 总 体 压 力 的 连接 条 件 


eth ~ d : ee 
定义 总 体 压力 : = p, + f Ê, dE ， 总 体 压 力 对 x 求 导数 有 : 
d d 5 d® és 
si Sd cat Oe | a a Ea 
dx p, dx dS, dx dx \dP,'% dé dx dx 


- dé ~p 曲线 变化 非常 平缓 ， 所 以 可 近似 采用 : 


er y d 1 ~ dp. 1 w 
EM o- f irae is 且 根 据 ee 


ox— f PER (25) 


dp, 5. " dé 
- i, d® dp, ~ dp 
(21) 可 得 : —(1-w)=-—£- c (26) 
H 可 得 (1-a) rs Ío T 
将 (26) RA (20) A: 
= d® | 
7 = 一 KU 一 +7 (27) 
V, í o) dx “Te 
HH (27) 整理 得 到 如 下 方程 : 
dm_ F, Yo (28) 
dx (I-a)KZ, (1-a) 
Xf (28) 在 [x,,% ERA RAS: 
Xp d® Xp v Vo 
afr e? _ _ (29) 
Xa dx Xa aie 5 
a 区 A x 区 7, 
® -®, = z : (30) 
A e i e 
和 ly 


oR (1-a) (oR (lo) 
由 (30) 或 (31), 4x, >m 9%, A: 
PP,=, (32) 

说 明 在 驱 蔡 前 缘 界 面 处 ， 有 总 体 压力 连续 。 


>< 4.3 驱 蔡 前 缘 界 面 的 压力 的 连接 条 件 
将 (30), (32) 整合 得 到 新 的 界面 条 件 : 


v =n + B 
©, =o, 


根据 (33) 揭示 的 总 体 渗流 速度 和 总 压力 在 驱 替 前 缘 界 面 的 连接 关系 ， 结 合 达 西 公式 (21) 可 以 
计算 出 气 和 水 的 渗流 速度 在 界面 处 的 连接 条 件 ， 根 据 p. = p,- p, 以 及 总 体 压力 定义 ， 可 计算 气 和 水 
压力 的 在 界面 处 的 连接 条 件 : 

(ve -GP ) f? = fon(w -GP?)+ fo 6. 
(vi +GP*) f? = fen(v +GP )+ FE FEB 


(33) 


dp‘ ~, dp? 

pit], ie eao paj peas oa 
~, ape ~, dp? 

piapia) RPE- pr pt) Beas 


4.4 驱 蔡 前 缘 两 侧 气 -水 两 相 流量 计算 
联 立 (22) 和 (31) 可 解 出 驱 替 前 缘 两 侧 的 总 的 渗流 速度 ，vw 。 


, (©, -©,) 1-0 | nK’? (1-0) 
0 em) (x, —x,) (二 二 元 本 
(1-0) Kit no)R (oo) KY n(1-@") KA? 
(五 一 a | B 可 Vo (x, — x) 
D-o, (1-0) ea)" (1-0) 
p- 


a n(x- x,) i Ce g y=) p Ca. 
(1-o")K*A" (1-@”)K’A? (I-o")K*A" (1-0) KA? 


将 (35) 代入 (21) 可 以 分 别 计算 出 内 veo 


dp: (Ss) ase +(p%-p’) gsin 


~a _ Jara 
View fo + 


dS“ Ox 
Sh zara “GP By) OSE pa GN . 
ead; a ) A +o pi) g sind 


K (36) RA (22) TUER, vto 


z dpi (S5 j 
v= jin AT 
b b 
p? = Fie? dp, (st) os’ +(p° p’)gsind 


rn dS? x 


) 5 讨论 


C1) 本 文 研究 结果 显示 在 驱 蔡 前 缘 界面 处 ， 流 体 压 力 和 渗流 速度 可 以 不 连续 ;流体 渗流 速度 不 


(37) 


(38) 


处 区 域 的， 两 相 体 很 容易 在 孔 喉 位 置 友 生 截断 ; 


(2) 对 于 突变 界面 处 流体 是 否 能 够 不 连续 分 布 ， 学 术 界 还 有 争议 ， 一 些 学 


件 ， 可 避免 采用 高 密 网 格 采 样 条 件 下 ， 得 到 满足 计算 结果 。 
(3) 大 多 数 情 况 下 ， 分 相 流 系 数 随 压 力 变化 较 小 ， 可 忽略 〈o=0 )， 


C4) 驱 蔡 形成 的 移动 界面 与 多 孔 介 质 物 性 突变 产生 的 界面 相 比较 ， 流 体 渗流 速度 和 压力 在 界面 
的 连接 条 件 特 征 不 同 ， 即 移动 的 突变 界面 的 连接 条 件 与 静止 的 突变 界面 的 连接 条 件 不 同 , 需要 注意 应 


L 连续 是 因为 渗流 速度 记录 的 是 流量 除 以 渗流 模 截 面积 ,即使 流体 分 子 的 真实 速度 在 界面 处 连续 , 但 渗 
- 流速 度 可 以 不 同 ， 当 流体 在 界面 处 连续 分 布 , 流体 压力 场 在 界面 处 应 连续 ,但 在 微观 上 ， 突 变 界面 所 


者 根据 传统 的 界面 条 
件 认为 流体 不 可 能 连续 分 布 。 我 们 认为 , WAPEN, 从 等 效 计 算 的 角度 , 流体 即使 在 突变 区 连续 分 布 ， 
但 分 布 范 围 很 狭窄 ,渗流 阻力 占 比 极 低 ， 采 用 本 文 提 出 的 能 够 反映 跨 间断 界面 渗流 规律 的 界面 连接 条 


用 条 件 ， 例 如 ， 本 文 所 建立 的 界面 条 件 不 适合 描述 跨 尺 度 ， 跨 类 型 介质 的 流量 的 计算 ; 


(5) 本 文 以 气 -水 两 相 为 例 ， 没 有 考虑 溶解 性 ， 如 果 考 虑 溶解 性 ， 因 物质 守恒 方程 的 变化 ， 突 变 


界面 条 件 也 随 之 变化 ， 本 文 只 考虑 了 一 维 情况 ， 对 于 多 维 界面 条 件 更 为 复杂 ; 
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